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Abstract 

A dilute Bose gas system is studied using the functional Schrodinger pic- 
ture theory. The ground state properties are obtained by solving the infinite 
dimensional Schrodinger equation variationally. It is shown that a shifted 
Gaussian trial wavefunctional enables us to calculate a higher order correc- 
tion, which corresponds to the fluctuation contribution from the condensate. 
The obtained term is compared with the quantum correction arising from the 
low energy 3 — > 3 scattering. 
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I. INTRODUCTION 

Recent successful demonstrations of Bose-Einstein condensatiorJHl have spurred a revival 
of interest in Bose gases. In order to gain a deeper understanding than the mean field result .@ 
it is necessary to calculate correction from quantum fluctuations around the mean field. The 
problem of interacting Bose gas has a long history. The standard result for the ground state 
energy density given in standard many-body textbookJFi has the form 



E„ 2nh 2 an 2 f 128 . „ 1/2 
— = < 1 H =(na 3 ) v + 
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where n = N/V is the number density of atoms and a is the S-wave scattering length. 
Although the coefficient of ln(na 3 ) in the na 3 correction was calculated early in 1950' 
a successful calculation of the constant k had to wait until recently. Braaten and Nieto0 
employed an effective field theory and the renormalization scheme to obtain the value of k in 
terms of a and atomic radius Rq. In this paper, we show that the interaction of particles out 
of or into the condensate give also an na 3 contribution to the ground state energy in addition 
to the 3 — > 3 scattering contribution of Braaten and Nieto.0 The relative importance of each 
term will be discussed. 

In the following, we apply the functional Schrodinger picture(FSP) scheme to a non- 
relativistic Bose gas to obtain the ground state properties. The Schrodinger picture for 
many particle systems or relativistic quantum field theory is an infinite dimensional extension 
of quantum mechanics and leads into a functional formulation.lil The FSP method in the 
relativistic quantum field theory is known to have several advantages over the conventional 
Green's function method.Eil'Ill One of such advantages is that the FSP method can readily 
be combined with a variational method to yield non-perturvative results. 

The FSP formulation of many-body fermionic system has been shown to yield the 
Hartree-Fock results when coupled with a Gaussian trial functional.ll!lll In order to obtain 
a result going beyond the mean field value, we employ a shifted Gaussian trial functional 
and calculate a contribution arising from the condensation fluctuation. 
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II. MANY-BODY DILUTE BOSE GAS 

We start with the standard model Hamiltonian for an interacting Bose gas 

H = Y ~2^T a fc afc + 2V ^ a I 1 a l 2 a fc3 a fe4^ 1 +fc 2 ,fc3+fc4) (2) 

k k i k 2 
k 3 k 4 

Here, the constant matrix element U should be determined by requiring that H produces 
the two-body scattering properties in vacuum. Since the system is a dilute gas, we assume 
that the scattering length a is much less than the inter-particle spacing n 1//3 . Therefore, 
the validity of the present calculation is restricted by the condition na 3 <C 1. Under this 
assumption, the zero momentum operators do, a\ become nearly classical due to occurrence 
of a condensate. The interacting part of the Hamiltonian is rewritten explicitly in terms of 
clq , as 

U t t 
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+ — ^2 \a\ +q a\a k a q + al +q al q a k a , (3) 

v k?&Q 

where only terms of order of A*"q , Nq, and y/No have been retained. Three zero momentum 
operators term does not exists, because any selection of three zero wave-vectors make the 
fourth zero by the momentum conservation. It should be noted that the last term arises 
from the interactions of particles out of and into the condensate, thus representing the 
condensation fluctuation.! 

In order to study the contribution from the last term, we introduce a new variable 7^ 

Ik = Yl a l+q a Q- ( 4 ) 

Also, we take the mean field approximation of 7^ as J2 q ^o( a k+q a q)- Replacing the operators 
a , Oq by V^O; H int approximately becomes 
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The number operator satisfies the relation 
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Using this relation and keeping up to the terms of v N, the model Hamiltonian becomes 
1 
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where e° = h 2 k 2 /2m. The ground state properties of the above Hamiltonian is to be calcu- 
lated using a variational scheme of the FSP. 



III. FUNCTIONAL SCRODINGER PICTURE METHOD OF BOSONIC 

SYSTEMS 

The FSP approach to many-boson system is simpler than many-fermion systemsBll 
when the Bose-Einstein condensation(BEC) is not included. However, when the BEC exists, 
a special care is necessary to include the condensation in the theory In the FSP theory, the 
bosonic field operators satisfy the equal time commutation relation 

*(x),¥(y)] =6{x-y). (8) 

Since we are working on a homogeneous system, it is convenient to work in the fc-space; 



We find that the fundamental commutation relations in Eq. (|8|) is satisfied, if a k and a\ are 
expressed as functional operators as follows, 



5<j)(k) 



(10) 



The weighting factor a k is necessary to correctly take account the BEC. a k is assumed 
real and even under under k — > —A;. The above expressions and the commutation relation, 
4>(k), —j^rj = S(k — k'), yield the constraint 

a k a k >5(k - k') = 5 kk >. (11) 



Note that a k 5(0) = 1 except k = 0. This is true only when a k is almost zero-flat for any k 
except k — 0. 

The Schrodinger equation in k space is expressed as 
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In the next section, we solve this functional differential equation using a variational method. 



IV. THE GROUND STATE ENERGY DENSITY 



The explicit form of the Hamiltonian in Eq. (|T2"| ) is obtained by substituting Eq. ( |TUD 
into Eq. (0), 
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(13) 



Since we are looking for the ground state, we may assume that the ground state functional 
is real, positive, and has no nodes anywhere.li! Thus, we write 



(14) 



where C is a normalization constant. When this expression is substituted into in Eqs. (|12| 
and (P^D , we obtain the relation 
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The power counting technique^ tells us that G[<p] is minimally quadratic in 4>(k). Thus, 
one can write G[4>] shifted Gaussian in d>, 

G[<f>] = J d 3 k [g{k)<f>{k)<t>{-k) + f(k)<f>(k)] , (16) 

where g(k) and f(k) are variational function parameters that will be determined in the 
below. The shifted Gaussian form of the trial function for the ground state allows us to 
calculate the contributions from linear terms of field operators. The trial functional method 
does not work if the Hamiltonian contains triplet field operator terms. This explains why 
we had to introduce the mean field operator 7^ in the Hamiltonian. 
Substituting Eq. (0) into Eq. (|i~5l), we obtain the following relation, 
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Counting powers of both sides, we obtain 
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and 



E g = \n 2 VU + [nUg(k) - (e° + nU) 



nU 

The constraint, Eq. ([TT|), was used above. 

Eq. ( |l8t ) for the variational function can be easily solved to yield that 

e° k + nU + E(k) 
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where 



E{k) = yJ{e° k + nU) 2 -{nU) 2 . 



(23) 



E(k) represents the excitation energy spectrum. 

The expressions of g(k) and f(k) allow us to evaluate the ground state energy. 
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The first two terms, Eq and £?i in the ground state energy are exactly same as those ob- 
tained through other methods in the literature. Following the standard step relating the 
scattering length to the interaction, 
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we obtain the ground state energy density of the first two terms as 
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The total ground state energy in the present calculation have an additional contribution 
E2, which will be calculated in the following sections. The corresponding wavefunctional is 
given by 



$ o [0] = CexplJ d 3 k 



2nU 
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7^ and a k are discussed in the following sections. 



V. THE PARTICLE DEPLETION 



(27) 



The particle depletion from the zero momentum condensate is defined by 

N-N 1 f 



N N 



(28) 



This value is equal to Jo/N in Eq. @, and can be calculated using Eqs. fllOD , fllil) , and 
(0). We obtain for =< a\a k >, 
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The average value of 4> 2 {k) and <p(k) in k state are also obtained as 
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Substituting the above results into Eq. ( p9|) , we obtain the average number of particle in 



state as 



n k = -1 + 



(32) 



where k = is not included. Clearly, we observe that a k cannot be 1 in the FSP formalism 
with the BEC. At the same time, n k should satisfy the following two conditions: 



limnfc = oo, and lim n k = 0. 

fe^O fe— *oo 



(33) 



We now use a simple comparison method to calculate the particle depletion. We have 
only two free variables of k, g(k) and E(k), in the ground state wavefunctional. f(k) does not 
contribute to the average particle depletion but represents the fluctuation of the condensate. 
From Eqs. ([H]) - (|23"D we obtain the following properties: 
(i) As k -> 0, 

1 



9(0) 



and E(0) = 0. 



(ii) As k — » oo, 



nUg(k) ~ E(k) (-eg). 



(34) 



(35) 



Using the above properties, we observe that the simplest way to satisfy these conditions is 

nUg(k) 



(36) 



2 E{k) ' 

Indeed, a k has the desired properties of a very sharp peak at k = 0, and almost zero values 
elsewhere. 

The particle depletion is now calculated as 
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The the summation over k and substitution of the S-wave scattering length a yield the 
standard resu 



N-N 



na 



3\ V2 



(38) 



N 3 V 7i J 

So far, we have shown that the FSP approach coupled with the variational method can 
produce the standard results on the ground state energy and the particle depletion in a 
straightforward manner. In the following, we discuss a higher energy correction. 
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VI. THE FLUCTUATION CORRECTION TO THE GROUND STATE ENERGY 

DENSITY 



The higher order contribution to the ground state energy in Eq. (E^) comes from the 
interactions of particles out of and into the condensate. Clearly, this term represents the 
fluctuation contribution from the condensate and disappears if the number of condensate 
particles are conserved strictly. The fluctuation contribution, E 2 , carries the ultraviolet 
divergence. The same divergence in Ex was handled through the expression of the scattering 
length a to order of U 2 as given in Eq. (p5|). However, for the calculation of E 2 such a simple 
cancellation scheme is not possible, because it is now necessary to expand a up to U 3 . 
Therefore, we resort to an altanative cutoff procedure, called "minimal subtraction" .0 In 
minimal subtraction, linear, quadratic, and other power ultraviolet divergences are removed 
as part of the regularization scheme by subtracting the appropriate power of k from the 
momentum space integrated. Following this prescription, we obtain E 2 in a dimensionless 
form 



E 2 2U 



k 



1 1 



e°, + 2nU -° 



(39) 



where E c = 2irh 2 an 2 V/m. 

Now, it is necessary to evaluate 'jkl-k = J2q,q'^o a l+q a q a -k+q' a g'- ^ i s known that a 
dominant contribution arises when a particle interacts with itself, which happens when 
q' = k + qJzl Therefore, 

lkl-k ^ a l+q a k+ q a\a q 

q^O 
q^O 



'na 3 



=n*-3jvv < 4o) 

In the second step, we used the fact that no is much larger than rik(k ^ 0). For the last 
step, Eqs. (^) and P5| ) are used. This approximation gives a simple result for the second 
order correction, E 2 , 
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This fluctuation correction has the same na 3 dependence with the 3 — » 3 scattering contri- 
bution obtained by Braaten and Nieto.0 This point will be discussed in detail later in this 
section. 

In the present approximation scheme, the total ground state energy density is given by 
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And the corresponding chemical potential is given by 
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It was shown that the low-density expansion for the energy density has the general 
form,i'0 
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(44) 



where [n/2] is the Gauss number which takes the integer part of n/2. Comparing this 
equation with Eq. (0), we find the standard results C 00 = 1, C w = 128/15y / 7r, and C21 = 
8(47r — 3\/3)/3. In addition, Braaten and Nieto calculated the constant under logarithm, 
which corresponds to Ki = 8 ( 4?r ~ 3v/ 3) ln(2.8i?o/ fl2 ); when the scattering length a is much larger 
than the atomic radius Rq. Since the result, ^2, we obtained within the approximation of a 
shifted Gaussian functional does not depend on the atomic parameters and has a different 
physical origin, we believe that the two constants are independent and additive so that 
k = Hi + k,2- Comparing the two terms, we observe that K2 is dominant if Ro/a < 0.75. In 
the real alkali atom gases, generally Ro/a ~ 0.1.0 Therefore, we believe that K2 is a larger 
contribution in general. 



11 



VII. SUMMARY 



We have presented a functional Schrodinger picture approach to a dilute Bose gas sys- 
tem. Employing a shifted Gaussian trial wavefunctional, we have calculated the ground 
state energy density and the particle depletion. It is shown that the present scheme allows 
to calculate the ground state properties in a straightforward manner and, in addition, to 
obtain the fluctuation contribution beyond the mean field result. The obtained correction 
is compared with the quantum correction term obtained through the low energy 3^3 
scattering. 

In the functional Schrodinger picture scheme, the many-body Schrodinger equation 
is solved directly without resorting to the Green's function or a canonical transforma- 
tion. Therefore, this scheme provides an intuitive alternative to conventional many-particle 
physics. 
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